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^^ ' Abstract. We calculate intensity and angular power spectrum of the cosmological 

*vj ' background of synchrotron emission from cold dark matter annihilations into electron 

positron pairs. We compare this background with intensity and anisotropy of 

astrophysical and cosmological radio backgrounds, such as from normal galaxies, radio- 

O ,. galaxies, galaxy cluster accretion shocks, the cosmic microwave background and with 

Galactic foregrounds. Under modest assumptions for the dark matter clustering we 
^ ' find that around 2 GHz average intensity and fluctuations of the radio background at 

^ , sub-degree scales allows to probe dark matter masses > 100 GeV and annihilation cross 

sections not far from the natural values (av) ~ 3 x 10~^^ cm^ s~^ required to reproduce 
the correct relic density of thermal dark matter. The angular power spectrum of the 
^t^ I signal from dark matter annihilation tends to be flatter than that from astrophysical 

Qs ' radio backgrounds. Furthermore, radio source counts have comparable constraining 
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1. Introduction 

Whereas it is known from cosmological observations that cold dark matter represents 
a fraction Vt^ ~ 0.233 of the total present energy density of our Universe [1, 2], its 
nature is still elusive. Dark matter can not only be detected directly in dedicated 
experiments searching for nuclear recoils from the scattering of dark matter particles, 
or produced in particle accelerators such as the LHC, but can also reveal its existence 
indirectly [3]: Although, apart from dilution from cosmic expansion, the density of 
dark matter does not change significantly after self-annihilations freeze out in the early 
Universe, residual self-annihilation can give rise to significant fluxes of 7— rays, electrons, 
positrons, neutrinos, and even some antimatter such as anti-protons and positrons, 
especially in regions with large dark matter densities. The energies of the secondary 
particles can reach up to the dark matter particle mass which can be of order a few 
hundred GeV. Secondary electrons and positrons can annihilate and give rise to a 511 
keV line emission, and they emit synchrotron radiation in the magnetic fields of galaxies 
which can be detected in the radio band. Therefore, cosmic and 7— ray detectors, 
neutrino telescopes, and even radio telescopes can be used for indirect dark matter 
detection as well. 

One of the most promising dark matter candidates are weakly interacting massive 
particles (WIMP) such as they are predicted within supersymmetric extensions of 
the Standard Model. Such particles have masses mx ^ 100 GeV. If they are 
produced thermally, in order to reproduce the correct average dark matter density, 
their annihilation cross sections have to be {av) ~ 3 x 10~^^ cm'^/s whereas large values 
are possible in case of non-thermal production. 

Traditionally, indirect dark matter detection has focused on high energy emission, 
specifically signatures in 7— rays [4, 5, 6, 7, 8, 9, 10, 11, 12, 13]. High energy 
emission of neutrinos, the particle the most difficult to detect, can be used to establish 
conservative constraints on the total annihilation cross section, which are of the order 
~ 10~^^cm^/s [14]. The fluxes of antiprotons and positrons from galactic dark matter 
annihilations has also been extensively used to constrain dark matter properties [15, 
16, 17]. Furthermore, the synchrotron radiation emitted by dark matter annihilation 
products close to the Galactic centre has also been studied. It has been found that if the 
dark matter proflle close to the central black hole is a spike formed by adiabatic accretion, 
typical dark matter annihilation cross sections within supersymmetric scenarios can lead 
to intensities comparable to the radio emission observed from the Galactic centre [18, 19]. 
It has further been shown that measurements of the radio flux away from the Galactic 
centre by the WMAP experiment strongly constrain the annihilation cross sections to 
values < 10~^^ cm^/s for mx — 100 GeV [20, 21]. A multi- wavelength analysis of dark 
matter annihilations from the Galactic centre has recently been performed in Ref. [22]. 

In the present paper we evaluate the diffuse synchrotron emission from the electrons 
and positrons produced by dark matter annihilation in the cosmological distribution of 
dark matter halos. We compute both its overall intensity and its angular power spectrum 
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as well as the distribution of visible dark matter annihilation sources as a function 
of apparent luminosity. We will find that comparing the resulting signals with other 
backgrounds and foregrounds under conservative assumptions allows to test annihilation 
cross sections close to the natural scale {av) ~ 3 x 10~^^cm^/s. 

In Sect. 2 we provide the general setup of our calculations. In Sect. 3 and 4, we 
apply it to astrophysical backgrounds and the dark matter induced signal, respectively. 
In Sect. 5 we compare the overall diffuse signal and its anisotropy to other foregrounds 
and discuss the resulting dark matter constraints, and in Sect. 6 we conclude. Finally, 
an appendix presents technical details of the calculations. We will use natural units in 
which c = 1 throughout. 

2. Setup 

We consider a distribution of sources which emit a radio luminosity per frequency 
interval L{h',V,z) which depends on a parameter V, on frequency u and on redshift 
z. The energy flux per frequency interval and solid angle is then given by 

^^'^'-J'^'dzdQj' dV^ ' ' 4iTdL{zy ' ^ ^ 

where for abbreviation we write Uz = {1 + z)v, {dn/dV)(V, z) is the co-moving volume 
density of objects per unit interval in the parameter "P, dL{z) is the luminosity distance, 
the factor 1 + z comes from redshifting the frequency interval dv, and the co-moving 
volume per solid angle and redshift interval is 
d'V ^ d,{zf ^ r_{zf_ 

dzdn (1 + zyH{z) H{z) ■ ^ ^ 

Here, for a flat cosmological geometry, the Hubble rate is 

H{z) = Ho [n„,(l + zf + ^a]'^' , (3) 

r{z) = /t(^) (1 + z)dt = Jdz'/H{z') is the co-moving distance and t{z) = Jq dz'/[{l + 
z')H{z')] is cosmic time as function of redshift. Throughout this paper we will assume 
a fiat, ACDM Universe with the total non-relativistic matter density Qnr = 0.279 and 
the dark energy density Q^ = 1 — Qnr — 0.701 (all other contributions to the energy 
density are negligible) in units of the critical energy density pc = 3HQ/{87iG]<i), where 
Gn is Newtons constant and Hq = H{0) = 100 h km s"^ Mpc"^ with h = 0.701 [1, 2]. 

In order to calculate the anisotropics we introduce an emissivity of squared power 
per frequency interval, £^(z/, k, z), 

C\v, k, z) = Cjiu, k, z) + Cl{v, k, z) , (4) 

which, similar to the approach in Ref. [9], we split into the two parts Cf{i',k,z) and 
jC,2{i^,k,z) and which also depends on the co- moving wavenumber k. The first part is 
essentially Poisson noise and corresponds to the sum over squared luminosities, 

Ci{,.,k,z) = I dV^{P,z)[L{,^z,V,z)\u{k,V)\]\ (5) 
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where in the following we define J^f{k) = J d^re^^''^ f (r) as the spatial Fourier transform 
of any function /(r) and where u{k, V) = J^u{k, V) is the Fourier transform of the spatial 
emission density u{r,V) of an individual source, normalized to unity, J d^ru{r,V) = 1. 
The second contribution to Eq. (4) is determined by the correlation between sources, 

1 2 



C2{iy,k,z) = Piin{k,z) 



J dV^{V,z)L{u,,V,z)b{r,z)\u{k,V)\ 



(6) 



where Piir^{k,z) = J d^re'^^'^{6p/p){r, z) is the linear power spectrum of the density 
fiuctuations {6p/p){r, z) and we have also introduced a bias factor b{V, z) of the sources 
with respect to the density field. 

The angular power spectrum Ci is given by 

Ci = {\aim\') , (7) 

where 

a,^ = J dn[Jiu,n) - {Jiu))]Y;^in) (8) 

in terms of the spherical harmonic functions Yi„i{Vl) and the intensity Jlu, Q) measured 
along direction Q. For a statistically isotropic sky this results in 

'^J dzda [iirddzn- ' ^ ' 

Using dii^z) = (1 + z)r{z) and inserting Eq. (2) in Eqs. (1) and (9) finally gives 

1 r dz r dn 

and 

[ATTy J dL{zYH{z) 

Formally, for point-like sources, the integral over redshift in Eq. (11) is divergent at 
z — > 0. In practice this is regularized by the fact that the nearest source has some 
minimal distance and that one can subtract the most luminous point sources which are 
also the nearest sources. In addition, the integral is regularized by the spatial extent of 
the sources, represented by the factor \u{k,V)\'^ in Eqs. (5) and (6). The role of these 
effects in practical calculations will be discussed in Sect. 5.2. 

3. Astrophysical Sources 

For astrophysical sources, V can be identified with the radio luminosity L^^ at some 
fixed frequency vq. Eq. (10) then simplifies to 

where Lcut(^) = AT:di,{z)^Scvx,l{^ + ^) is the intrinsic luminosity corresponding to the 
apparent point source fiux ^cut above which we consider the source to be resolvable and 
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thus subtractable from the diffuse background. For the multipoles Eqs. (5) and (6) can 
then be written as 

n2 



dL,, 



C^{u,k,z) = / dL^^L 

and 

icut(2) 



' ^"^ iK,z) 



Liu,) 



u{k, z) 



Cl{u,k,z 



P\m{k, z) 



"-^I/Q-^fO IT 



dn .^ \L{vz) ,, ,,,^ , 
[L^^, z)— — u{k, z)b{L,,^, z) 



(13) 



(14) 



respectively. For the luminosity functions dn/dLy^^ of normal and radio galaxies we will 
use the expressions given in Ref. [23]. 



4. Dark Matter Annihilation 

For annihilation of dark matter with mass mx and phase space averaged annihilation 
cross section times velocity (af), V can be identified with the mass M of dark matter 
halos. We then follow the approach of Ref. [18] and write 

{av) 



L{u, M) 



2m\ 



S{u,M) 



(15) 



where we define Eiv^M) as a quantity which does not depend on annihilation cross 
section or mass of the dark matter particles, 

o3 , ,ra^ Y^{> E) 



S{iy,M) 



VSe- 



nip 



d'rpi{v)B{v) 



mx 



dE- 



V 



Vc{E) 



(16) 



Psyn(^)+Pic(^)' 

In Eq. (16), e and vrif^ are the electron charge and mass, respectively, Ph(r) is the dark 
matter halo density profile, Biv) is the local magnetic field strength, , Fe(> -S) is 
the multiplicity per annihilation of electrons and positrons with energies larger than 
E, Psya{E) = 2e^-B^-E^/(3mg) = (IGe^n / 3)ubE'^ / m^ is the total synchrotron emission 
power of one electron of energy E in a, magnetic field of strength B, corresponding to an 
energy density ub = -B^/(87r), and Pic{E) = (16e'^7r/3)M^£'^/mg is the energy loss rate 
at energy E due to inverse Compton scattering on a low energy photon field of energy 
density u^. Furthermore, we use the function 

(17) 



F{x) = X K5/3{y)dy, 

J X 

in Eq. (16), with the critical frequeny 

Ait nie KnieJ 
In the following we use the approximation [24] 

F{x) ~ S[x - 0.29] 

such that Eq. (16) can be simplified to 



9 



£{iy,M) ~- 



m": 



1/2 



Ye[> E,{u)] 



8 V0.297r 



,1/2 



/(M) 



(19) 



(20) 
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where 

and the critical energy Ec{v) is the inversion of Eq. (18), 

In Eq. (20) we neglect the magnetic field dependence of Fe[> E(,{y)\. For mx ^ 100 GeV, 
B > a. few micro-Gauss and i^ ~ 1 GHz, the parameters we are interested in, this is a 
good approximation because the critical energy Eci^z) ^ mx/10. Typical values for 
these parameters are Yf. ~ 10 [18]. This corresponds to a fraction /e ~ 0.3 of the total 
annihilation energy going into pairs. The energy fraction going into pairs of energy 
above E can be expressed in terms of Ye{> E) as 

fe{E) = ^ J^"" dE'E'^iE') < 1 . (23) 

With the above expressions we can rewrite Eq. (10) as 
j(j,) = M J_ f^:!±_\ f d^ [dM^{M,z)Y,[> E,(z/,)]/(M).(24) 

Furthermore, we can redefine C\ and C\ from Eqs. (5) and (6) by extracting constant 
factors and write 

C\{u, Kz) = j dM^{M, z) (F,[> E,{v,)]I{M)\u{k, M)\f (25) 

and 

Cl{v, k, z) = PUK z) (j dM-^{M, z)Y,[> E,{u,)]I{M)b{M, z)\u{k, M)\\ , (26) 

where u{k,M) relates to the halo profile, u{k, M) = J dhe^'pl{r){eB)-^/'^{r)/I{M). 
With these quantities we can now write 

' 1024-0.297r3z/ l^m2,y 7 {1 + z)dL{zyH{z) ' ^ ' 

Details about the quantities that enter these expressions are given in Appendix A. 

Eq. (27) can also be obtained as follows: Limber's equation relates the two- 
dimensional angular power spectrum ^2(0 ^o ^^e three-dimensional power spectrum 
P^{k) in the fiat sky approximation [25]: Given a three-dimensional statistically random 
field /(r) = f{Q,r), one considers the observation at r = of the projection 

P{n)= f drw{r) f{n,r) (28) 

with some given radial weight function w{r), where r is the co-moving distance. If the 
field / fiuctuates on scales much smaller than the characteristic scale over which w{r) 
varies, then we have 

Qc^ jdr'^Pjil/r.zir)) (29) 



Dark Matter Signatures in the Anisotropic Radio Sky 



where Pf{l/r,z) is the power spectrum of {f{Qi,r)f{Q2,r)) at the co-moving 
wavenumber k = l/r. 

Neglecting the variation of the magnetic field B within the halo regions contributing 
most to the annihilations, the radio intensity Eq. (24) along a given direction Vt can be 
written as 



J{iJ,n) 



{av) 9p^ 



m': 



1/2 



,Al±^y^^, E,.J1^3^ , (30) 



m\ QAttu^/^ \0.29TTeB J J ^" H{z) ^^'^ "^"^'^ 1 + u^/ub 

where pm = ^mPc is the average dark matter density at zero redshift, and 6 = Sp/p 
is the relative overdensity. Because the dominant contribution comes from the dark 
matter halos, where 6 ^ 1, we can approximate (1 + 5)^ ^ 6"^. Assuming a constant B 
and a constant optical photon field of density u^p — 5eVcm~'^, we can write the factor 
(1 + u^/ub)~^ = [1 + ""op/mb + Mo(l + zY /ub] , where Uq is the CMB energy density 
at 2; = 0. This factor effectively cuts off the redshift integration at 2; ^ 2. Since Ec{vz) 
varies little over this redshift range, we can then further simplify Eq. (30) to 



J{u,n) 



Fe[> E,{v)] {av) 9p 



m" 



1/2 



m^ 647r\/eS^ V0.297r/ 

Comparing this with Eq. (28), we can use 

for the random field and the weight function is 

Fe[> E,{iy)] (av) 9pi 



dz- 



{l + zf'H\z,n) 



Hiz) 



l + ^ + ^(l + z)4 



w{z) 



■^m 



m" 



1/2 



:i + zf/^ 



777,1 



0.297r 



^ + ^(1 + Z)4 
UB UB^ ' 



.(31) 



(32) 



.(33) 



(■x 647rvefii 

The power spectrum P/(A;, z) appearing in Eq. (29) is then the Fourier transform of the 
two-point correlation function of / in real space. Following Ref. [8], Pp.(k^z) can be 
written as the sum of a one-halo and a two-halo term, Ps^ik, z) 
with 



plhl 

-^52 



TDlht 



nMcut(2) 



{k,z)= / dM—[J^s^{k,M,z)f= / 



Psf{k,z)=P,Uk) 



PUk) 



dM 

Mcutiz) 

-'■*'-'- min 

M^ut(z) 

IVI Yn, i n 



dM 



dn 
dM 

dn 



r-Mcut(z) ^^dn 
iM„,„ dM 

1 2 



\k,z)+Pl^{k,z 
'AbTAk,M,z)\" 



piri^+^r 



(34) 



dM-—b(M) 
dM ^ ' 



b{M)Ts^{k,M,z) 

AhJ^Ak,M,z) 



P^(l + ^) 



(35) 



where Mmin is the minimal halo mass and Mcut{z) is the halo mass corresponding to the 
apparent point source flux Scut above which we consider the source to be resolvable and 
thus subtractable from the diffuse background. Furthermore, Afy is a boost factor which 
accounts for possible substructure in the halos. The average of the clumping factor 
appearing in Eq. (31) is given by 

At 



6\z) 



Pm(l + Z) 



*'-^""'(^) ,.. dn 

dM-^_ X 



M„ 



dM 



dV,,pi{T,M,z) 



(36) 
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where dVh is the halo volume element. 

A generic form for the halo mass function dn/dM appearing in the equations above 
was first proposed by Press & Schechter (PS) [26]; a modified version of this form is 
given by Sheth and Tormen [27] (ST). When comparing the results obtained from these 
two forms, we find differences by factors less than 2. Thus, we adopt the PS formula 
throughout our paper. 

Current knowledge of the dark matter density distribution mostly comes from N- 
body simulations, and the universal dark matter profile firstly proposed is the Navarro- 
Frenk- White (NFW) model [28]. Combining Eq. (A. 14) with Eq. (A. 17) in Appendix A, 
in this model the dark matter profile within each halo can be written as 

^ Ac(^) C^ Pm(^) .07N 

^''^''^ 3 ln(l + c)-c/(l + c) r/r,(l + r/r,)2' ^ > 

where r^ is a characteristic radius, and the concentration parameter c is defined as the 
ratio of r^ and the virial radius r^, c = r^/rg, see also Appendix A for more details. 
Note that r^ is not a free parameter, but depends on M and c because r„ is related 
to M via M = 4TTAc{z)pm{z)/3, where Ac{z) in an Einstein-de Sitter Universe is 
about 187r^. With the above definition, c and M completely determine the dark matter 
distribution of a given halo. The minimal halo mass is still rather uncertain. The value 
^min = 10"^ Mq [29] is close to the free-streaming mass [30, 31, 32], below which there 
are no fluctuations in the dark matter density to form a halo. Note that the magnetic 
field may be much smaller than micro Gauss scales in such small halos. In contrast, the 
value Mjnin = 10^ Mq roughly corresponds to the minimal mass of dwarf galaxies which 
are known to contain micro Gauss scale magnetic fields [33]. We, therefore, choose 
^min = 10^ Mq as fiducial value in the following, noting that the dark matter signal 
would increase by only a factor about two for Mmjn = 10~^ Mq. We will furthermore 
use B = 10/iG as fiducial value for the magnetic field. This is a realistic value given that 
most annihilations occur in the densest regions where also magnetic fields are somewhat 
larger than typical average galactic fields. 

The clumping factor is very sensitive to the concentration parameter, namely oc c^. 
N-body simulations indicate that the concentration has a log-normal distribution [34] 
with a median value of 

c{M, z) = 4 i±^ , (38) 

1 + z 

where the collapse redshift z^ is implicitly given by the relation M^{zc) = O.OIM, where 
M^[z) is the mass scale at which o"(M^,,z) = 6c- How the concentration parameter 
depends on halo mass and redshift is still an open question. One can extrapolate Eq. (38) 
to minimal halo masses M ~ 10^^ Mq. When comparing the parameterization Eq. (38) 
with high resolution simulations [29] we find that it gives realistic values for the minimum 
halo mass. It is a conservative estimate because at z ~ it gives values c ~ 70 for the 
minimum halo mass which is significantly smaller than other parameterizations [35, 36]. 
Recent studies show that dark matter halos exhibit considerable substructure [37, 
38, 39, 40, 41, 42, 43, 44, 45, 46, 47]. The total mass of these substructures only account 
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for about 10% of the host halo, but they can give an extra boost factor Ai, ~ 10 for dark 
matter annihilation. Some studies show that if one takes into account substructure and 
assumes a cuspy center slope [48, 49, 50], the theoretical prediction can well explain 
the excess of high energy positrons and the diffuse 7— ray background observed by 
the Heat [51, 52] and EGRET [53, 54] experiments, respectively. The subhalos follow 
a certain mass and redshift distribution which is still unknown. Therefore, to be 
conservative we assume the NEW halo model and simply parametrize any possible boost 
factor with the parameter A^ ~ 10. The substructures occur on small scales and do not 
influence the power spectrum in the range we are interested, / < 10^. 

5. Results 

5.1. Diffuse Radio Emission 

The cosmic microwave background (CMB) dominates the radio sky at frequencies above 
~ 1 GHz, whereas astrophysical sources such as normal galaxies and radio galaxies 
dominate at lower frequencies down to kHz frequencies [23]. Recently it was argued that 
synchrotron emission of strong intergalactic shocks can also significantly contribute to 
the diffuse extragalactic radio below 500 MHz [55, 56]. 

Using the formuli developed in Sect. 4, we now evaluate the contribution of 
synchrotron emission from pairs produced by dark matter annihilation in the magnetic 
fields of dark matter halos. We consider neutralinos as dark matter candidate, and for 
the following figures we assume a neutralino mass of 100 GeV and a total annihilation 
cross section of {av) = 3 x 10~^^ cm^/s, fixed for reproducing the correct relic density for 
thermal relics. We also assume that the average total number of electrons and positrons 
per annihilation is Yg — 10, and that the halo substructure implies a boost factor 
Af) ~ 10. We compare the resulting dark matter signal with astrophysical contributions 
to the diffuse background that can be computed from the expressions in Sect. 3. 

Eor astrophysical sources the diffuse radio background is likely dominated by normal 
galaxies and radio galaxies. To estimate the contributions from these sources, we follow 
Ref. [23], which use the observed correlation between the radio and infra-red flux of 
galaxies. This approach assumes that the radio emission is related to the star formation 
and is sensitive to the redshift evolution of the sources, but can explain the observed 
radio background quite well. 

EoUowing the above assumptions, in Eig. 1, we show the different contributions to 
the average diffuse radio intensity. Eor astrophysics sources, normal galaxies contribute 
more than radio galaxies. This is because although the individual radio galaxy is brighter 
than a normal galaxy on average, this is overcompensated by the larger number of normal 
galaxies. Also shown in Eig. 1 is a possible contribution from intergalactic shocks [55, 56] 
normalized such that its angular power spectrum is comparable to the one of the Galactic 
foreground, see Sect. 5. 2. 

Of course, the CMB absolutely dominates the radio sky in the wide range from v ~ 
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Figure 1. The average diffuse background flux intensity with no point-source removal. 
Contributions from normal galaxies (blue curve), radio galaxies (red curve), fr-om radio 
and normal galaxies combined (black curve), and from a scenario for radio emission 
from galaxy cluster shocks (magenta curve) [56] (see text for the normalization) are 
compared to our fiducial dark matter annihilation scenario with mx = 100 GeV, 
{av) - 3 X IQ-"^^ cu)? /s, At = 10, B = 10 /iG, !/„„„ = 10^ AIq (brown curves). 
Here, the solid brown curve is for Y^. = 10, while the dashed brown curve is for 
Ye{E) ~ nix/E. Also shown is the CMB background (cyan solid curve) as well as 
its subtractable part, determined by uncertainties of the absolute CMB temperature 
(dotted cyan curve). The Galactic foreground at Galactic latitude b > 20° is shown as 
the green band within uncertainties. 



1 GHz to a few hundred GHz [57], and above these frequencies Galactic foregrounds such 
as dust emission dominates. Since the CMB is a black body radiator its contribution to 
the solid angle averaged radio flux can be subtracted up to the uncertainty of its average 
absolute temperature. Currently the CMB temperature is measured to 2.725 ± 0.001 
K [58]. We convert this temperature uncertainty into an intensity of CMB confusion 
noise. Fig. 1 shows that this confusion noise dominates other astrophysical backgrounds 
and the diffuse signal of our fiducial dark matter scenario at z/ > 4 GHz. At lower 
frequencies the dark matter signal vJ{v) tends to decrease as ^/v for Y^ ~ const., see 
Eq. (31), whereas the background from normal galaxies tends to be flat, see Fig. 1. There 
is thus an optimal window at frequencies i^ ~ 1 GHz where dark matter annihilation 
signatures can be detected and where self-absorption is negligible. Constraints on dark 
matter parameters can, therefore, only be established for annihilation cross sections 
about a factor ten higher than the fiducial cross section required for thermal dark matter. 
In addition, there are three diffuse foregrounds from our Galaxy in the frequency 
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range we are interested: The first is synchrotron radiation emitted by high energy 
electrons gyrating in the Galactic magnetic field, the second is free-free emission from 
the thermal bremsstrahlung from hot (> lO'^K) electrons produced in the interstellar 
gas by the Galactic UV radiation field, and the third foreground is dust emission which 
arises from the thermal re-radiation of absorbed stellar light. Fig. 1 shows that these 
foregrounds tend to dominate the astrophysical backgrounds and the dark matter signal 
in the fiducial scenario. 

Can we test the properties of dark matter more powerfully? The absolute CMB 
temperature is difficult to measure more precisely than to the current permille level, 
because of inevitable systematic errors. Small-scale temperature fiuctuations AT/T ~ 
10~^ have been seen by the COBE and WMAP satellites because temperature differences 
can be measured more precisely since systematic errors cancel in measurements of 
temperature differences. Furthermore, if the Galactic foregrounds have a smooth 
directional dependence, they may pose less of a contamination when considering the 
anisotropy of the radio sky. We, therefore, consider in the following the angular power 
spectra of the radio sky in order to see if it can provide further tests of dark matter 
properties. 

5.2. Anisotropy 

Whereas the diffuse average radio flux provides only one number at a given frequency 
to compare with other astrophysical and cosmological backgrounds, potentially much 
more information is contained in the angular power spectrum. For example, the power 
spectrum as a function of angular scale tends to be different for dark matter annihilation 
and astrophysical sources because the contribution of the latter to the diffuse radio flux 
is dominated by fewer bright sources. Our goal in this section is whether this can provide 
dark matter signatures or constraints on mass and annihilation cross section. 

Before calculating the angular power spectra, we discuss their qualitative behaviors. 
The angular power spectrum Ci = C}'^ + Cf'^ can be divided into one-halo {C}^) and 
two-halo {Cf^) terms, corresponding to the two contributions Eqs. (5) and (6) to Eq. (4), 
and thus to Eq. (11). The two-halo term arises from the correlation between distinct 
halos which is described by the linear power spectrum. The one-halo term represents 
correlation within the same halo. Both one-halo and two-halo term are proportional 
to !«(/;;, P)!^, the square of the Fourier transform of the spatial emission proflle. At 
large angular scales, \u{k,V)\ ~ 1, such that C}^ is essentially independent of /. The 
one-halo term is thus sometimes called Poisson noise. At scales comparable to the 
size of the source, \u{k,V)\'^ starts to become suppressed. Therefore, both one-halo 
and two-halo terms are expected to be suppressed for multipoles / larger than the 
typical distance to the source divided by the linear source size. The two-halo term 
is furthermore proportional to the linear power spectrum which is also suppressed for 
co-moving wavenumbers k > 0.03 Mpc^^. Therefore, the ratio of the two-halo term to 
the one-halo term is suppressed for / > 0.03Mpc~^rH — 100, where th — 3000 Mpc is 
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the Hubble scale. The one-halo term eventually dominates at very small angular scales. 

In Eq. (11), for point-like sources, formally the one-halo term C}^ would diverge 
for 2;min — > 0, whereas the two-halo term Eq. (6) is regularized by the linear power 
spectrum Piin(/c, 2), which is suppressed at large k = l/r{z). This is because the 
flux of nearby sources of a given luminosity diverges. We can ignore such sources 
because they can be identified as individual bright sources and be removed from the 
background flux in actual observations. We can remove sources with intrinsic luminosity 
Lc\it{z) > At: d L^zyS cut/ {)- + z), corresponding to the point-source sensitivity S'cut of 
the telescope. Alternatively, one can regularize Eq. (11) by integrating from some 
finite minimum distance corresponding to the typical distance to the nearest source, 
'^min ~ 1 Mpc. Furthermore, Eq. (11) is also formally regularized at z^^^ — > by the 
spatial extent of nearby sources, described by \u{k^V)\'^. For the NEW profile, the 
mass of the halo within distance r from the halo centre increases as r^ up to r = r^, 
and then increases logarithmically between r^ and r^ since Ph{f) oc r""^, see Eq. (37). 
Therefore, the dominant contribution to the halo mass comes from r < rg. Similarly, 
for r < rg the annihilation signal increases as r, but between r^ and r„ increases only 
as r~^ — r~^. Assuming the emission traces ph for astrophysical emission processes 
and p\ for dark matter annihilation, the Fourier transforms of these dependencies 
then give u{k,V) oc k~'^ for k ^ r~^, with 7 = 2 for astrophysical emission and 
7 = 1 for dark matter, see Appendix A. 6 for more details. Since k = l/r{z), and 
thus \u [l/r{z),V] p oc r(z)^''', the one-halo term in Eq. (11) diverges only for 7 < 0.5. 
Therefore, under our assumptions for the emission profile, Eq. (11) is convergent even 
without cut-offs in either r^^i^ or the apparent luminosity. Since nevertheless in particular 
the one-halo term is quite sensitive to nearby sources, in the following we study its 
dependence on S'cut and r^i^. 

Fig. 2 shows the dependence of the two-halo term on S'cut- According to Eq. (6), the 
two-halo term scales with the square of the average flux. Since the apparent luminosity of 
radio galaxies can be of the order of a Jansky (1 Jy = 10~^^ erg cm~^ Hz~^ s~^), the two- 
halo term from radio galaxies starts to decline when we cut sources being less luminous 
than a critical luminosity below a Jansky. In contrast, the contribution of normal 
galaxies which are much less luminous than radio galaxies starts to decline only when 
we cut sources more luminous than ~ lO/xJy. The contributions of dark matter halos to 
the dark matter annihilation signal is basically unaffected by any source removal, even 
to luminosities down to ~ Ip Jy. This is easy to explain: In our fiducial scenario the 
largest dark halos of about 10^^ Mq produce only about 1.3 x 10^*erg/s at 2 GHz from 
dark matter annihilation, far less than the typical radio luminosity of galaxies of about 
2 X 10^°erg/s. As a result, removing bright sources increases the contribution of dark 
matter annihilation to the two-halo term relative to the contribution from astrophysical 
sources. 

Next we discuss the one-halo term. The one-halo term is more sensitive to the 
cut-offs in apparent luminosity S'cut and to the minimal distance rmin than the two-halo 
term because of two reasons: First, the two-halo term Eq. (6) is the square of an integral 
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than Scut to the two-halo term at 2 GHz. The red, blue and brown lines represent the 
contribution from radio galaxies, normal galaxies, and dark matter (fiducial scenario 
with Ye = 10), respectively. 



of luminosities, whereas the one-halo term Eq. (5) is essentially Poisson noise and thus 
proportional to an integral of squared luminosities, which makes the contribution from 
bright sources more important. Second, the two-halo term is further regularized by the 
linear power spectrum at large k = l/r{z). In Fig. 3 we show the cumulative contribution 
of sources dimmer than Scut to C/'*. Similarly to the two-halo term shown in Fig. 2, the 
contribution of radio galaxies and ordinary galaxies decreases rapidly below ^ 1 Jy and 
10/iJy, respectively, whereas the contribution of dark matter annihilation is affected 
less by source removal. Nevertheless, the contribution of bright sources is now much 
larger than for the two-halo term, as expected, and the one-halo term continues to rise 
with inclusion of brighter sources. On the other hand, practically one should cut off the 
integral at some minimal distance r^^^ ~ 1 Mpc within which there are essentially no 
bright sources. Since sources at small distance appear bright, the cut-off in luminosity 
and minimal distance is of course to some extent degenerate, as confirmed by Fig. 3. For 
radio galaxies, removal above ^ 10^ Jy is equivalent to restricting to distances larger 
than 1 Mpc. For ordinary galaxies, cutting at a minimal distance rmin = 1 Mpc is 
equivalent to removing sources brighter than 0.1 Jy. Since observational sensitivities 
are considerably better than these luminosities, cutting at r^im — 1 Mpc does, therefore, 
not introduce any significant uncertainties. Note that in Fig. 2 dark matter dominates 
the two-halo terms if all sources above ~ 0.1/i Jy are removed, while in Fig. 3 it would 
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Figure 3. The cumulative contriloution of sources of apparent luminosity S smaller 
than S'cut to the one-halo (Poisson) term at 2 GHz. The solid and dotted curve 
represent the cases of Tmin = and Tmin = 1 Mpc, respectively. Color keys are as 
in Fig. 2. 



dominate the one-halo terms only for unrealistically small cut-off luminosities < 1 nJy. 
This is because the one-halo term is much more sensitive to bright sources than the 
two-halo term and because the dark matter contribution consists of dimmer sources 
than ordinary astrophysical sources. 

The angular power spectra of the radio background at 2 GHz produced by galaxies 
and by our fiducial dark matter scenario are shown in Figs. 4, 5, and 6 for different source 
removal cuts. Based on the above discussion, the qualitative behavior of the one- and 
two-halo terms can be easily understood: In Fig. 4 we assume galaxies to appear point- 
like and we remove sources brighter than 0.1 mJy. The one-halo terms from these sources 

1 /9 

thus increase proportional to [/(/ + 1)] in the above figures. The same apphes to the 
one-halo term of the dark matter contribution for / < 10^, corresponding to angular 
scales 9 ^ vr// > 0.02°. At smaller angular scales the power spectrum is suppressed 
by the inner structure of the dark matter halos. We can estimate this critical scale as 
follows: The one-halo term is dominated by the brightest halos which correspond to the 
largest and nearest halos. In our fiducial scenario, the annihilating dark matter in the 
largest halos can emit a radio flux of ~ 7 x lO^^ergs"^ Hz~^ at 2 GHz, such that the 
minimum co-moving distance is r ~ 830 Mpc, z ~ 0.2 if sources brighter than 0.1 mJy 
are removed. The scale r^ for the corresponding 10^^ M© halo is about 0.21 Mpc. This 
corresponds to a multipole / ~ vrr/r^ ~ 1.2 x 10^. This simple estimation is consistent 
with our detailed calculation shown in Fig. 4. 
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Figs. 4, 5, and 6 show that for radio galaxies the one-halo term is always larger 
than the two-halo term at all multipoles, as expected because of the high luminosity 
of radio galaxies. For dark matter and normal galaxies, the two-halo term dominates 
at small /. The dependence of the angular power spectrum on / can potentially be 
used to discriminate the dark matter signal from astrophysical contributions: For 
/ < 3 X 10^, the annihilation power spectrum looks significantly flatter than the signal 
from normal galaxies. In other words, at large angular scales, the annihilation signal has 
relatively more power. This can be understood as follows: After cutting bright sources, 
many more dim nearby annihilation sources than galaxies contribute. In addition, at 
large redshift the synchrotron emission from dark matter annihilation is suppressed by 
the increased inverse Compton scattering rate on the CMB, see Eq. (33). The two- 
halo term is proportional to Pwr^ik) which peaks at ~ 0.03 Mpc~^, corresponding to 
/ ~ 0.03r(z)/Mpc. The on average smaller distance to the dark matter halos then 
translates into relatively more power at small /. 

In Fig. 5, we take into account the spatial extent of the radio emission of galaxies. 
We assume the luminosity profiles of galaxies to be roughly proportional to the dark 
matter density profile which is obtained following Appendix A. 2 with the halo mass is 
obtained from the relation between mass and bolometric luminosity [59]. As a result, for 
normal and radio galaxies the one-halo term starts to drop for / > 6000 and / > 2.5 x 10"^, 
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Figure 7. Angular power spectra of the radio sk;y at 2 GHz compared with various 
estimates of tlie Galactic foreground at Galactic latitude h > 20° (green shaded region) 
and the CMB (cyan curve). The brown band represents the annihilation spectrum, 
where the upper and lower ends correspond to -Fdm = 10 and Fdm = 1, respectively, 
see Eq. (39), and from which halos brighter than 0.1 mJy were removed. The black- 
dotted and black-solid curves represent the total signal from normal and radiogalaxies, 
for luminosity cuts S'cut = lOmJy and Scut = 0.1 mJy, respectively. Also shown is a 
possible contribution from intergalactic shocks [56], normalized such that its angular 
power spectrum is comparable to the Galactic foreground. 



respectively. Compared to the dark matter signal, the suppression thus sets in at slightly 
smaller / for normal galaxies, but only at larger / for radio galaxies. For normal galaxies 
this is due to the more extended emission profile which more closely follows the density 
as opposed to the squared density in case of dark matter. This is also reflected by 
the Fourier transform of the emission profiles shown in Fig. Al. For radio galaxies this 
effect is overcompensated by the fact that they are much brighter such that after cutting 
bright nearby sources, their average distance is much larger where their angular extent 
appears smaller. 

Since future radio detectors such as the square kilometer array (SKA) [60] can 
reach point flux sensitivities of ~ 1/xJy, we show the power spectra of the background 
remaining after a corresponding luminosity cut in Fig. 6. Since the two-halo term from 
dark matter annihilation is insensitive to such luminosity cuts whereas the contribution 
from galaxies decreases rapidly, as shown in Fig. 2, the relative contribution of dark 
matter annihilation increases and gives rise to a flatter power spectrum at moderate /. 

We now have to compare the cosmological background power spectra discussed so 
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far with other potential contaminations. Fig. 7 compares the signals from ordinary and 
radio galaxies and from our fiducial dark matter scenario with the power spectra of the 
CMB and of the Galactic foreground at high Galactic latitude. The power spectrum 
of the Galactic foreground is not very well measured and we represent its uncertainties 
as a green band in Fig. 7. At high Galactic latitude below 10 GHz local foreground 
fluctuations dominate over the CMB power spectrum which is why the CMB anisotropy 
measurements are performed above 20 GHz. Concerning annihilation signatures of dark 
matter with mass mx ^ 100 GeV in the angular power spectrum of the radio sky, 
the optimal frequency band is around 2 GHz. At higher frequencies, the synchrotron 
emission of electrons produced from dark matter annihilations cuts off due to Eq. (22) 
and the CMB signal increases. At lower frequencies, synchrotron emission by Galactic 
electrons dominates the power spectrum even at high Galactic latitude [61, 62, 63, 64]. 
Around 2 GHz, Galactic synchrotron emission always dominates, whereas free-free 
emission is a factor few smaller. 

Also shown in Fig. 7 is a possible signal from intergalactic shocks [56]. Since its 
normalization is rather uncertain, we normalized it such that it is comparable to the 
average estimate of the Galactic foreground. The thermal SZ effect [65] is another 
characteristic contamination caused by hot ionized gas in galaxy clusters and filaments 
outside of clusters [66]. Since it dominates at small angular scales, / > 3000, and at 
high frequencies above 30 GHz, we can neglect this effect here. 

As can be seen from comparing Fig. 5 and 6 and from Fig. 7, future radio telescope 
arrays sensitive around z/ ~ 2 GHz, with their higher point flux sensitivities should 
allow to further reduce the contribution from galaxies, whereas for / < 6000 the dark 
matter contribution is hardly changed by removing still fainter sources. This can be 
understood from the fact that the dark matter signal is dominated by the two-halo term 
which is insensitive to Scut for Sent ^ 1 nJy, see Fig. 2. This shows that for dark matter 
annihilation the distribution of /(/ + l)Ci is nearly fiat for 200 < / < 2000. At smaller / 
the power spectrum is dominated by Galactic foregrounds and at larger / the one-halo 
term from galaxies grows rapidly. The most sensitive range 200 < / < 3000 should be 
accessible to present and future radio telescopes with their high angular resolution. The 
SKA will have a sensitivity of about 6 x 10~^^ ergcm"^ sr~^ s~^ in the units of the above 
figures. 

We conclude that the power spectrum from dark matter annihilation tends to be 
flatter than other contributions because of an interplay of the following effects: 

• The astrophysical signals are dominated by fewer and much brighter sources than 
the dark matter annihilation signal which consists of many faint sources. For 
'S'cut = 0.1 mJy, the two-halo term dominates for / < 10^ for both the signals 
from galaxies and from dark matter annihilation. In addition, the dark matter 
signal is significantly flatter in that angular range, i.e. it has relatively more power 
at small /. This is because the two-halo term is proportional to the linear power 
spectrum whose peak in wavenumber for the on average closer and dimmer dark 
matter annihilation sources translates into smaller / at these luminosities. 
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• For / > 10^, the inner spatial structure of the galaxies and dark matter halos 
becomes important. The inner structure tends to suppress the power spectra, but 
the exact angular scale at which these effects become important depends on the 
halo size, the profile of the emission and source luminosity cut-off. 

• The various components evolve differently with the Universe expansion. For 
example, at high redshift inverse Compton scattering on the CMB tends to suppress 
synchrotron emission in dark matter halos, whereas astrophysical sources such as 
radio galaxies tend to be more active at z ^ 3. 

5.3. Dark Matter Constraints 

We can now scale the dark matter signal to parameter values different from the fiducial 
scenario, by multiplying with the factor 

'M (Yl\ I <^^) ^ /lOOGeVy xio^x 1/2 

lOy UOy Isx 10-26cm3s-V V mx ) \ B ) ^ ' 



dm 



We caution that a boost factor as high as A^ ^ 10 has not been verified in all dark matter 
structure simulations and that the average magnetic field B could be significantly smaller 
than 10/iG if many small-scale subhalos contribute. However, smaller values for A}, and 
B partially compensate in Eq. (39) so that one could still obtain observable signatures 
as long as B is large enough to produce emission at GHz frequencies, see Eq. (22). 

If we choose Fdm = 10 (upper end of brown band in Fig. 7) and Scut = 0.1 mJy, the 
annihilation spectrum dominates over other cosmological backgrounds for 100 < / < 10^ 
and should become distinguishable from the Galactic foreground. Note that this 
foreground is likely further reduced close to the Galactic poles. In this situation it 
should thus be possible to disentangle the rather flat power spectrum of the dark matter 
annihilation signal from other contributions in the range of 200 < / < 3000. We can 
thus assert that radio observations are sensitive to 

i^dm>10, (40) 

with some dependence on the source luminosity cut-off S'cut- Note that the dark matter 
signal shown in Fig. 7 does not strictly scale with Fdm because it depends on sources 
dimmer than S'cut? here chosen as 0.1 mJy. However, since few brighter dark matter 
halos contribute, the signal scales with Fdm in first approximation. 

One can also compare observed radio source counts as a function of apparent point 
source flux with predictions for astrophysical sources and dark matter annihilation 
sources. This is done in Fig. 8 for the same parameters as used in Fig. 7. This establishes 
the constraint Fdm < 10. In contrast. Fig. 7 provides dark matter signatures for future 
measurements but currently does not allow to put a constraint on Fdm because of the 
uncertainties in the Galactic foreground spectrum. Note that a future observational 
extension of the source count spectrum in Fig. 8 to apparent luminosities S < /iJy 
will provide an additional test for dark matter which predicts a shallower source count 
distribution than astrophysical sources. 
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Figure 8. Observed radio source counts {dN/dS)S'^'^ as function of apparent radio 
flux S compared with predictions for normal galaxies (blue curve), radio galaxies (red 
curve), and annihilations from dark matter halos (brown band, for 1 < i^dm ^ 10). 
Green shaded region and triangles are data from Ref. [67] . 



We have not computed the contribution from dark matter annihilations in our 
own Galaxy to the anisotropic radio flux in the present work. However, we know from 
Refs. [10, 21] that for our fiducial values for cross section and mass, at least the smooth 
halo component does not lead to fluxes higher than current observations from WMAP. 
The contribution from Galactic substructures is probably more model dependent than 
our cosmological flux which apart from an overall boost factor depends only on the host 
halo distribution and effectively averages over a much larger ensemble of halos. This 
can also be seen from Ref. [68] where the predictions of the 7— ray flux from Galactic 
dark matter annihilations varied over orders of magnitude. 



6. Conclusions 

Many different indirect detection signatures have been investigated to constrain the 
parameter space for dark matter [8, 12, 69, 70, 71, 72, 73]. In the present paper we have 
calculated intensity and angular power spectrum of the cosmological background of 
synchrotron emission from the electrons and positrons produced in annihilations of cold 
dark matter. The resulting radio background around ~ 2 GHz and its angular power 
spectrum for multipoles 200 < / < 3000 has comparable or better sensitivity to dark 
matter annihilation cross sections than other signatures. Furthermore, a comparison 
of observed radio source counts with predictions for dark matter annihilation results 
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in the constraint -Fd^ < 10 for the parameter defined in Eq. (39). Under reasonable 
assumptions on dark matter clustering and magnetic fields in the halo environment, 
the range of annihilation cross sections corresponding to the constraint Fdm ^ 10 
is comparable to constraints from synchrotron emission in an NFW profile [10, 21]. 
Galactic 7— ray constraints derived under similar assumptions [74] are also comparable. 
The sensitivity of our signal is considerably better than conservative limits based on 
annihilation into neutrinos [14], and comparable to limits on annihilation into 7— rays 
from diffuse cosmological emission [13]. 

Sensitivities to values of order ten for the parameter F^^ defined in Eq. (39) 
are interesting for non-thermal dark matter whose annihilation cross sections can be 
larger than our fiducial value {av) = 3 x 10~^^cm^s~^, the cross section required 
for thermal dark matter. An example for a non-thermal dark matter candidate with 
large cross section is the wino LSP occurring in supersymmetric theories with anomaly 
mediation [75, 76]. We believe that radio observations in particular with future 
instruments such as SKA can provide valuable information on dark matter. 
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Appendix A. Cosmological Dark Matter Distribution 

Appendix A.l. Power Spectrum and Halo Mass Function 

The mass function of the halo distribution is derived from the Press-Schechter 
formalism [26]. In this approach fluctuations in the linear density field with 5 > 5c 
decouple from the local Hubble expansion of the universe and collapse to form non- 
linear structures. The fraction of the volume that has collapsed is predicted to be 

/eou(M(i?),^) = '^ rdSe-'^/"^'^^'^^ , (A.l) 

V27rcr(it, Z) J 5c 

where R is the co-moving radius over which the density field has been smoothed, which 
is related to the halo mass by M{R) = pm4:7iR^/3 with pm the co-moving matter density 
of the universe. The number density of halos is then found to be given by [26, 27] 

dM ~ M dM ~M^^^'dM' ^ '' 

where 



/(^) = V ^^[1 + K')-1e-- . (A.3) 
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Here 

and 

5c{z) - 1.686 (A.5) 

is the critical density required for spherical collapse at a redshift z in an Einstein-de 
Sitter space. The variance in the density field smoothed with a top-hat filter of radius 
R= (3M/47rp^)^ is 

a\M,z) = G\z)j'^^^^^^\W{kR)\\ (A.6) 

where 

3 

W{x) = ^-[sinfx) — xcos(x)] , (A. 7) 

Piinik) is the hnear matter power spectrum, and 

H{z)jrdz'{i+z')[H{z')r^ 

^^'^ - H,J-dz'{l + z')mz')]-^ ^^-'^ 

is the growth factor with hnear perturbation theory, often also denoted hj D{z). In 
Eq. (A. 3) A, p, and a are constants, with the canonical Press-Schechter (PS) and 
Sheth-Tormen (ST) mass functions corresponding to the parameters {p = 0,a = 1) 
and {p = 0.3, a = 0.707), respectively. The normalization A is determined by requiring 

mass conservation such that 

1 /""^ dn r°° 

— dMM—= duf{u) = l. (A.9) 

Pm Jo aM Jo 

For PS A = 1 and for ST A = 0.3222. 

The primordial power spectrum P{k) oc A^k""' can be modified by the content and 

evolution of different matter components of the Universe due to the perturbations that 

enter the horizon at different epochs. This allows one to relate the linear power spectrum 

to the primordial power spectrum through a transfer function T{k) via 

PiUk, z) = D\z)PUk. z = Q) = D\z)A,{k ■ M^cY^T\k) (A.IO) 

Fitting formula for an adiabatic CDM model give [77] 

TcDM(g) = ^''^V34g^^^^ ^^+^-^^^+^^^-^^^'+^^-^^^^'+^^-^^^^'^"'^'^^-^^^ 

where q = k ■ Mpc/ (hT) and F = Q,mh exp[Q,b{l + J{2h)/Vlm)\- One usually uses the 
rms fluctuation on an 8/i~^ Mpc scale to normalize the amplitude of the present power 
spectrum. From WMAP 5 year data, we adopt ris = 0.96 [1, 2], and Ag = 1.4 x 10^. 
Following Ref. [78], we may furthermore write the linear growth factor as 

D{z) = ^—^, (A.12) 

^^ 1 + ^^(0)' ^ ^ 

where an approximate expressionfor g{z) is 

^^ ' nrn{z)^/'^-nA{z) + {i + nUz)/2){i + nA{z)/70)- ^- ' 
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The halo mass function has to be supplemented by the dark matter density profile. For 
the dark matter profile within each halo we use an NFW profile [28], 



Phir) 



(A. 14) 



r/rs{l + r/rsy 

where r^ is a characteristic radius. Within the context of the spherical collapse model, 
the outer extent of the cluster is taken to be the virial radius 



Tv 



3M 



47rp^(l + zf^c{z) 



(A.15) 



where Pm(l + zY is the average physical background matter density of the universe at 
redshift z, and 



^ciz) 



ISvr' 



1-a 



215 \n^(l + zf^ 



(A.16) 



is the overdensity of the halo relative to the background density [79]. The ratio of the 
virial radius to the scale radius is called the concentration parameter c = r^jr^. A 
combination of the definitions of virial mass and density profile gives 

^^'^" 3 ln(l + c)-c/(l + c) ^^-^^^ 

Together, c and M completely determine the dark matter distribution of a given halo. 
Note that an NFW profile is a conservative assumption compared to steeper profiles 
that have been proposed, for example the Moore profile [80], which would consequently 
lead to a larger dark matter signal. 



Appendix A. 3. Concentration Distribution 

It is still uncertain how the concentration parameter depends on mass and redshift. In 
our paper we consider a realistic value of c{M = 10~^ Mq,z = 0) ~ 70. Extrapolating 
to the low mass range, we use [34] 

c{M,z)=a\±^, (A.18) 

1 + z 

where the collapse redshift z^ is implicitly given by the relation M^{zc) = O.OIM, where 

M*(z) is the mass scale at which a{M^:, z) = 6c- In a less conservative parametrization 

motivated by numerical simulations c would follow a log-normal distribution with 

standard deviation ac = 0.18, 



V{\nc\M,z) 



exp 



[In c - In c(M,z)]^ 



where the mean concentration parameter c is related to the halo mass via [36] 



c{M, z) 



Co 

l + z 



M 
M,{z = 0)_ 



(A.19) 



(A.20) 
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where Cq and a^ are constants whose numerical values [34] are typically chosen to be 
Co = 9 and Uc = 0.13. However, application of this parameterization to low-mass 
halos and to high redshift give values inconsistent with some simulations [29]. A third 
parameterization [35] uses 

ibiz) 



c{M, z) = a{z) 



M 



(A.21) 



M,{z) 

with a{z) = 10.3(l + z)"°-^ , and 6(^) = 0.24(l+z)-°-3. The parametrization Eqs. (A.19), 
(A. 20) give comparable dark matter signals, whereas the parametrization Eqs. (A.21) 
would lead to signals about a factor 20 higher than our conservative calculation. 

Appendix A. 4- Bias 

For the linear dark matter halo bias b{M, z) appearing in Eq. (6) we adopt [81] 

MM,,^l.^, (A.22) 

whereas for the galaxy bias we simply use unity. 

Appendix A. 5. Bolometric Luminosity 

In order to compute the Fourier transform of galaxy density profiles, we first need the 
galaxy mass. For normal galaxies we use the relation [59] 

M = 14^^^^>^Mo, (A.23) 

where Lgo is the luminosity at 60 microns, and z/go ~ 5000 GHz is its frequency, and 
Lq ~ 3.9 X lO^^erg/s is the solar bolometric luminosity. For radio galaxies we simply 
adopt their typical mass of about 10^^ Mq to estimate the Fourier transform of the their 
density profile. 

Appendix A. 6. Fourier Transforms 

The Fourier transform of the spherically symmetric NEW profile of mass M can be 
written as 

^,,(A;,M) = £p,{M,rf-^^^47rr'dr, (A.24) 

and analogously for JF,2(A;, M). For the purpose of plotting these Fourier 
transforms, see Fig. Al, it is convenient to renormalize them to unity for k — *• 
by introducing the new functions yi{k,M) = J^p^{k,M)/M and y2{k,M) = 
J^p2{k,M)/ JdVhplir). We then have yi{0,M) = 1, and yi{k > 0,M) < 1 for 
i = 1,2. For the NEW density profile, J dVhplir) = fcMpm^dz), where fc = 
(cV9) [1 - (1 + c)-3] / [log(l + c) - c/(l + c)]^and A{z) ~ 200. 

For the NEW profile, the mass of the halo within radius r increases oc r^ for r < rg, 
and then increase logarithmically for r^ < r < r^ where Ph{f^) oc r^^. Therefore, the 
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Figure Al. The normalized Fourier transforms yi{k,M) (solid lines) and y2(k,M) 
(dotted lines) of ph and p/j2, respectively, as functions of co-moving wavenumber k. 
The vertical lines denote the scale k = l/rs(M). 



dominant contribution to the halo mass comes from r < rg. Similarly, the annihilation 



signal is produced mainly within r 



< 



increasing there oc r, but increases only 



oc rl 



for r<r< 



Fig. Al shows that for kr^ ^ 1 we have yi^2 
y2{k, M) oc k-^, and yi{k, M) oc k~^ . 



1, whereas for kr^ ^ 1 one has 



Appendix A. 7. Foregrounds 

The radio intensity !„ at a given frequency v can be expressed in terms of antenna 
temperature T^(i^) via I^ = 2z/^/cb^a(i^)/co, where cq is the speed of light. Alternatively, 
Ii, can be written in terms of the thermodynamic temperature as the temperature 
of a blackbody with the given intensity at frequency u, thus Ii, = 2u^/{e^ — 1), 
where x = hu/k^T with h the Planck constant. Thus, for power law spectra 
/^ oc i^°, T4 oc i/"~^. In general, the CMB is expressed in terms of thermodynamic 
temperature T, while Galactic and extragalactic foregrounds are expressed in term of 
antenna temperature. Thermodynamic and antenna temperature are then related by 
T = TA(e^ - l)/x, and their fluctuations by AT = ATA(e^ - lf/{x^e''). For the CMB, 
X = hv/{k^TcMB) ^ z//(56.8GHz) with the CMB temperature Tcmb = 2.725 K [58]. 
Since we consider frequencies z/ < 10 GHz in the present paper, x <^1 and thus T c^Ta 
and AT ~ AT^. 
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From the definition of T^ we get 



-25 



-2s-^Hz~^sr-^ 



(A.25) 



/. = 3.06 X 10-1 ^J ^-J erg cm 

Since for Ii, oc i^" the power spectrum Ci" of I^ at frequency v scales as z/^", we can 
express it in terms of the power spectrum C^^{v') of the antenna temperature T^ at 
frequency v' via 



C/-'(z/) = 3.06x 10 



-25 



a-2 



^^s~^Hz"Sr 



1 „„-l 



(A.26) 



Here, a = —0.9 for synchrotron emission and —0.15 for free-free emissions, 
respectively [62]. For the normalization and the dependence on /, we adopted the 
best-fit model from observations at 2.3 GHz [63]. These parametrizations have been 
used in Fig. 7. 
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Figure A2. Diffuse energy spectrum of 7— rays from dark matter annihilation for 
our fiducial scenario, with Mmin = 10~^ Mq. This is consistent with Fig. 1 in Ref. [8] 
within about 10%. 



Appendix A. 8. Diffuse Flux and Power Spectrum of'-f—rays 

We can use our approach also to compute diffuse energy spectra and angular power 
spectra of 7— rays from dark matter annihilation. For neutralinos the spectrum of 
7— rays of energy E^ per annihilation can be parameterized by the simple expression [4] 



dN^ 

He 



{E, 



0.73 



g-7.776B^/mx 



mx (^^/mx)^-^ + 0.00014 



(A.27) 
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Figure A3. Angular power spectrum of 10 GeV 7— rays from dark matter annihilation 
in our fiducial scenario. Brown lines are for minimal halo mass M,„i,i = 10^ Mq, 
whereas red lines are for M„iin = 10^^ Mq. The higher and lower curves at high I 
denote the one-halo and two-halo terms, respectively. Note that C; is now normalized 
to the total intensity squared, I^. This is consistent with Fig. 6 in Ref. [8] within a 
factor ~ 2. 



To reasonably match the EGRET data [53], and compare with the results in Ref. [8], we 
multiply the predicted average 7— ray intensity spectrum by the boost factor Ai, ~ 240 
due to substructure within the host halos. Note that this boost factor is much more 
extreme than A;, ~ 10 assumed in the present work. We then have 



w{E^, z) 



(av) fQ 



dN^{E^,z) 



{l + zfE., 



(A.28) 



Svr \mx J ^ ' ■' ' dE^ 

for the weight function in Eq. (28). For our fiducial dark matter scenario we then obtain 
Fig. A2 for the solid angle averaged energy spectrum and Fig. A3 for the angular power 
spectrum at 7— ray energy E^ = 10 GeV. Note that if fainter sources were included, 
corresponding to smaller minimal halo mass Mmin, the relative fluctuations normalized 
to /^, and especially the one-halo Poisson term, would decrease. 

The diffuse energy spectrum is consistent with the results in Ref. [8] within about 
10%. In addition, the angular power spectrum at 7— ray energy 10 GeV is also consistent 
with the results in Ref. [8] within a factor two. This serves as an important cross-check 
of our results with independent calculations. 



Dark Matter Signatures in the Anisotropic Radio Sky 28 

References 

[1] E. Komatsu et al. [WMAP Collaboration], arXiv:0803.0547 [astro-ph]. 

[2] J. Dunkley et al. [WMAP Collaboration], arXiv:0803.0586 [astro-ph]. 

[3] G. Bertone, D. Hooper and J. Silk, Phys. Kept. 405, 279 (2005) [arXiv:hep-ph/G4G4175]. 

[4] L. Bergstrom, J. Edsjo and P. UUio, annihilations," Phys. Rev. Lett. 87, 251301 (2001) 

[arXiv:astro-ph/0105048]. 
[5] P. UUio, L. Bergstrom, J. Edsjo and C. G. Lacey, look," Phys. Rev. D 66, 123502 (2002) 

[arXiv:astro-ph/0207125]. 
[6] J. E. Taylor and J. Silk, Mon. Not. Roy. Astron. Soc. 339, 505 (2003) [arXiv:astro-ph/0207299]. 
[7] D. Elsaesser and K. Mannheim, background," Phys. Rev. Lett. 94, 171302 (2005) [arXiv:astro- 

ph/0405235]. 
[8] S. Ando and E. Komatsu, Phys. Rev. D 73, 023521 (2006) [arXiv:astro-ph/0512217]. 
[9] S. Ando, E. Komatsu, T. Narumoto and T. Totani, Phys. Rev. D 75, 063519 (2007) [arXiv:astro- 

ph/0612467]. 
[10] D. Hooper, G. Zaharijas, D. P. Finkbeincr and G. Dobler, WMAP Haze," arXiv:0709.3114 [astro- 
ph]. 
[11] A. Cuoco, J. Brandbyge, S. Hannestad, T. Haugboelle and G. Miele, Phys. Rev. D 77, 123518 

(2008) [arXiv:0710.4136 [astro-ph]]. 
[12] S. Ando, E. Komatsu, T. Narumoto and T. Totani, Mon. Not. Roy. Astron. Soc. 376, 1635 (2007) 

[arXiv:astro-ph/0610155]. 
[13] G. D. Mack, T. D. Jacques, J. F. Beacom, N. F. BeU and H. Yuksel, arXiv:0803.0157 [astro-ph]. 
[14] H. Yuksel, S. Horiuchi, J. F. Beacom and S. Ando, Phys. Rev. D 76, 123506 (2007) [arXiv:0707.0196 

[astro-ph]]. 
[15] L. Bergstrom, J. Edsjo, M. Gustafsson and P. Salati, JCAP 0605, 006 (2006) [arXiv:astro- 

ph/0602632]. 
[16] T. Bringmann and P. Salati, Phys. Rev. D 75, 083006 (2007) [arXiv:astro-ph/0612514]. 
[17] T. Delahaye, R. Lineros, F. Donato, N. Fornengo and P. Salati, Phys. Rev. D 77, 063527 (2008) 

[arXiv:0712.2312 [astro-ph]]. 
[18] G. Bertone, G. Sigl and J. Silk, Mon. Not. Roy Astron. Soc. 326, 799 (2001) [arXiv:astro- 

ph/0101134]. 
[19] G. Bertone, G. Sigl and J. Silk, Mon. Not. Roy. Astron. Soc. 337, 98 (2002) [arXiv:astro- 

ph/0203488]. 
[20] D. Hooper, Phys. Rev. D 77, 123523 (2008) [arXiv:0801.4378 [hep-ph]]. 
[21] P. Grajck, G. Kane, D. J. Phalen, A. Pierce and S. Watson, arXiv:0807.1508 [hep-ph]. 
[22] M. Regis and P. UUio, arXiv:0802.0234 [hep-ph]. 
[23] R. J. Protheroe and P. L. Bicrmann, Astropart. Phys. 6, 45 (1996) [Erratum-ibid. 7, 181 (1997)] 

[arXiv:astro-ph/9605119]. 
[24] Rybicki G.B., Lightman A. P., Radiative Processes in Astrophysics, 1979, John Wiley & Sons. 
[25] Kaiser, N., Astrophys. J. 388, 272 (1992) 
[26] W. H. Press and P. Schechter, Astrophys. J. 187, 425 (1974). 
[27] R. K. Sheth and G. Tormen, Mon. Not. Roy. Astron. Soc. 308, 119 (1999) [arXiv:astro- 

ph/9901122]. 
[28] J. F. Navarro, C. S. Frcnk and S. D. M. White, Astrophys. J. 462, 563 (1996) [arXiv:astro- 

ph/9508025]; 
[29] J. Diemand, B. Moore and J. Stadel, Nature 433 (2005) 389 [arXiv:astro-ph/0501589]. 
[30] S. Hofmann, D. J. Schwarz and H. Stocckcr, Phys. Rev. D 64, 083507 (2001). 
[31] X. 1. Chen, M. Kamionkowski and X. m. Zhang, Phys. Rev. D 64, 021302 (2001). 
[32] A. M. Green, S. Hofmann and D. J. Schwarz, JCAP 0508, 003 (2005). 
[33] P. P. Kronbcrg, Rcpt. Prog. Phys. 57, 325 (1994). 
[34] J. S. Bullock et al., Mon. Not. Roy. Astron. Soc. 321, 559 (2001) [arXiv:astro-ph/9908159]. 



Dark Matter Signatures in the Anisotropic Radio Sky 29 

[35] U. Seljak, Mon. Not. Roy. Astron. Soc. 318, 203 (2000) [arXiv:astro-ph/GGG1493]. 

[36] A. Cooray and R. K. Shcth, Phys. Rcpt. 372, 1 (2002) [arXiv:astro-ph/0206508]. 

[37] G. Tormcn, A. Diafcrio and D. Sycr, Mon. Not. Roy. Astron. Soc. 299, 728 (1998). 

[38] A. A. Klypin, S. Gottlobcr and A. V. Kravtsov, Astrophys. J. 516, 530 (1999). 

[39] B. Moore, S. Ghigna, F. Govcrnato, G. Lake, T. Quinn, J. Stadel and P. Tozzi, Astrophys. J. 524 

(1999) L19. 
[40] S. Ghigna, B. Moore, F. Governato, G. Lake, T. Quinn and J. Stadel, Astrophys. J. 544, 616 

(2000). 
[41] V. Springel, S. D. M. White, G. Tormen and G. Kauffmann, Mon. Not. Roy. Astron. Soc. 328, 

726 (2001). 
[42] A. R. Zentner and J. S. BuUock, Astrophys. J. 598, 49 (2003). 
[43] G. De Lucia et al, Mon. Not. Roy. Astron. Soc. 348, 333 (2004). 
[44] A. Helmi, S. D. M. White and V. Springel, Phys. Rev. D 66, 063502 (2002). 
[45] L. Gao, S. D. M. White, A. Jenkins, F. Stoehr and V. Springel, Mon. Not. Roy. Astron. Soc. 355 

(2004) 819. 
[46] L. Shaw, J. Weller, J. P. Ostriker and P. Bode, Astrophys. J. 646, 815 (2006). 
[47] J. Diemand, M. Kuhlen, P. Madau, M. Zemp, B. Moore, D. Potter and J. Stadel, arXiv:0805.1244 

[astro-ph] . 
[48] X. J. Bi, Nucl. Phys. B 741, 83 (2006) [arXiv:astro-ph/0510714]. 
[49] Q. Yuan and X. J. Bi, JCAP 0705, 001 (2007) [arXiv:astro-ph/0611872]. 
[50] J. Lavalle, Q. Yuan, D. Maurin and X. J. Bi, arXiv:0709.3634 [astro-ph]. 
[51] S. W. Barwick et al. [HEAT Collaboration], Astrophys. J. 482, L191 (1997) [arXiv:astro- 

ph/9703192]. 
[52] S. Coutu et al., Prepared for 27th International Cosmic Ray Conference (ICRC 2001), Hamburg, 

Cermany, 7-15 Aug 2001 
[53] A. W. Strong, L V. Moskalenko and O. Reimer, Astrophys. J. 613, 956 (2004) [arXiv:astro- 

ph/0405441]. 
[54] P. Sreekumar et al. [EGRET Collaboration], Astrophys. J. 494, 523 (1998) [arXiv:astro- 

ph/9709257]. 
[55] E. Waxman and A. Loeb, Astrophys. J. 545, Lll (2000) [arXiv:astro-ph/0007049]. 
[56] U. Keshet, E. Waxman and A. Loeb, Astrophys. J. 617, 281 (2004) [arXiv:astro-ph/0402320]. 
[57] G. F. Smoot, arXiv:astro-ph/9902201. 
[58] J. C. Mather, D. J. Fixsen, R. A. Shafer, C. Mosier and D. T. Wilkinson, Astrophys. J. 512, 511 

(1999) [arXiv:astro-ph/9810373]. 
[59] L. Spinoglio, M. A. Malkan, B. Rush, L. Carrasco and E. Recillas-Cruz, Astrophys. J. 453, 616 

(1995) [arXiv:astro-ph/9506139]. 
[60] C. A. Jackson, New Astron. Rev. 48, 1187 (2004) [arXiv:astro-ph/0409180]. 
[61] M. Tegmark and G. Efstathiou, arXiv:astro-ph/9507009. 
[62] M. Tegmark, D. J. Eisenstein, W. Hu and A. de Oliveira-Costa, Astrophys. J. 530, 133 (2000) 

[arXiv:astro-ph/9905257]. 
[63] G. Giardino, A. J. Banday, P. Fosalba, K. M. Gorski, J. L. Jonas, W. O'MuUane and J. Tauber, 

arXiv:astro-ph/0103233. 
[64] L. La Porta, C. Burigana, W. Reich and P. Reich, arXiv:0801.0547 [Unknown]. 
[65] R. A. Sunyaev and Y. B. Zeldovich, Astrophys. Space Sci. 7, 3 (1970). 
[66] Persi. F. M., Cen, R., Ostriker, J. P., Astrophys. J. 442, 1 (1995) 
[67] J. J. Condon, Ap. J. 338, 13 (1989). 
[68] J. M. Siegal-Gaskins, arXiv:0807.1328 [astro-ph]. 

[69] X. L. Chen and M. Kamionkowski, Phys. Rev. D 70, 043502 (2004) [arXiv:astro-ph/0310473]. 
[70] L. Zhang, X. L. Chen, Y. A. Lei and Z. G. Si, Phys. Rev. D 74, 103519 (2006) [arXiv:astro- 

ph/0603425]. 
[71] L. Zhang, X. Chen, M. Kamionkowski, Z. g. Si and Z. Zheng, Phys. Rev. D 76, 061301 (2007) 



Dark Matter Signatures in the Anisotropic Radio Sky 30 

[arXiv:0704.2444 [astro-ph]]. 
[72] Z. Myers and A. Nusser, arXiv:0710.0135 [astro-ph]. 

[73] N. Padmanabhan and D. P. Finkbcincr, Phys. Rev. D 72, 023508 (2005) [arXiv:astro-ph/0503486]. 
[74] S. Dodclson, D. Hooper and P. D. Serpico, Phys. Rev. D 77, 063512 (2008) [arXiv:071 1.4621 

[astro-ph]]. 
[75] L. Randall and R. Sundrum, Nucl. Phys. B 557, 79 (1999) [arXiv:hep-th/9810155]. 
[76] G. F. Giudice, M. A. Luty, H. Murayama and R. Rattazzi, JHEP 9812, 027 (1998) [arXivihep- 

ph/9810442]. 
[77] J. M. Bardeen, J. R. Bond, N. Kaiser and A. S. Szalay, Astrophys. J. 304 (1986) 15. 
[78] S. M. CarroU, W. H. Press and E. L. Turner, Ann. Rev. Astron. Astrophys. 30 (1992) 499. 
[79] J. P. Henry, Astrophys. J. 534 (2000) 565 [arXiv:astro-ph/0002365]. 
[80] B. Moore, T. Quinn, F. Governato, J. Stadel and G. Lake, Mon. Not. Roy. Astron. Soc. 310, 1147 

(1999) [arXiv:astro-ph/9903164]. 
[81] H. J. Mo, Y. P. Jmg and S. D. M. White, Mon. Not. Roy. Astron. Soc. 282, 1096 (1996) 

[arXiv:astro-ph/9602052]. 



